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A method  has  been  developed  for  performing  accurate  total y^iergy  calculations 
using  the  X-alpha  theory  without  the  usual  arbitrarinesa/in  selection  of 
calculational  parameters.  The  use  of  relatively  high/^values  and  atomic 
sphere  radii  in  the  ratio  of  their  atomic  (covalent)  radii  led  to  results  in 
agreement  with  experiment  for  dissociation  energies  and  equilibrium  separation 
distances  in  NO,  CO,  and  SN.  In  addition,  a method  of  obtaining  dipole  moments 
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I . INTRODUCTKM 


There  is  at  present  considerable  interest  in  the  causes  and  effects 

of  erosion  of  gun  tubes^.  Recent  experiaental  results  on  the  interaction 
of  combustion  gases  with  steel  alloys  indicate  that  chemical  effects 
can  be  important  in  the  erosion  of  steel  and  may  even  exceed  those  of  the 

2 

melt  and  wipe-off  mechanism  for  certain  combustion  conditions  . 

A reliable  means  of  cosqmtationally  simulating  erosion  processes 
would  be  a valuable  tool  to  aid  in  interpreting  the  experimental  data 
and  in  providing  estlMtes  for  erosion  behavior  under  conditions 
difficult  or  impossible  to  simulate  experimentally.  In  addition,  an 
accurate  microscopic  model  would  advance  our  fundamental  understanding 
of  hot  gas-metal  surface  reaction  mechanisms. 

A detailed  microscopic  or  quantum-chemical  treatment  of  the 
interaction  of  a gas  with  a metal  surface  may  be  divided  into  two 
separate  areas:  1)  a description  of  the  movement  of  the  atoms  during 
collisions  in  response  to  the  forces  exerted  upon  them;  2)  a 
description  of  the  forces  acting  upon  the  atoms.  This  latter  descrip- 
tion is  equivalent  to  specifying  the  potential  energy  of  the  system 
since  the  force  is  always  derivable  from  a potential  for  any  system  in 
which  energy  is  conserved. 

The  interaction  of  a gas  with  a metal  surface  is  intimately  involved 
in  the  process  of  chemisorption,  a process  which  may  be  loosely  defined  as 
the  adsorption  of  atoms  or  molecules  iqwn  a surface  with  a binding  energy 

in  excess  of  a certain  energy,  say  one  eV.^  In  many  cases  the  adsorption 
energies  are  high  enough  that  the  process  corresponds  to  electron  sharing, 

1. e.,  chemical  bonding. 

In  view  of  the  fact  that  computational  methods  of  quantum  chemistry 
are  having  increasingly  far  ranging  applications,  a task  was  initiated 
with  the  objective  of  determining  if  these  methods  might  be  profitably 
applied  to  understanding  chemisorption,  thus  enabling  one  to  gain 
insight  into  erosion  processes. 

During  the  execution  of  this  task  much  previously  unreported 
information  regarding  the  quantua-chemical  ni^el  was  obtained  and  is 
reported  here. 

TI  fH-Serwice  Oun  Yt*e  Veer  and  Erosion  Symposium,  US  Any  AnMment 
Research  and  Developswnt  Comnand,  Dover,  HJ,  Kerch  1977. 

2.  A.C.  Alkidas,  L.H.  Caveny,  M.  SuesMrfield,  and  J.W.  Johnson,  High 
Pressure  and  High  Teq>erature  Gas  Metal  Interactions,  presented  at 
ISth  JANNAF  Conbustimi  Conference,  Monterey,  CA,  September,  1976. 

3.  S.K.  Lyo  and  R.  Corner,  Interactions  on  Metal  Surfaces.  Springer- 
Verlag,  Berlin,  1975,  p.  4l. 
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II.  THEORY 


In  order  to  st^ply  useful  infonstion  on  the  process  of  cheni- 
sorption,  s theoretiesl  nodel  should  be  cspsble  of:  1)  describing 
excited  stntes,  2)  predicting  rensonnble  dissociation  energies, 

S)  predicting  reasonable  equlllbriuB  Intemuclear  separations:  and, 
for  soae  probleas,  4)  predicting  accurate  dipole  nonents.  Also, 
the  Mthod  should  be  citable  of  reasonably  priced  calculations.  An 
accurate  description  of  excited  states  is  necessary  since  nany  of 
the  atoas  will  be  in  excited  states.  The  need  for  an  accurate 
dissociation  energy  and  equlllbrluB  Intemuclear  separation  can  best 
be  understood  with  the  aid  of  Figure  I.  Two  atoas,  foming  a aolecule 
in  its  ground  state,  will  have  their  nuclei  separated  by  an  average 
distance  r^,  at  which  the  electronic  energy  will  be  a ninlaua.  Any 

separation,  r,  less  than  r^  will  cause  the  energy  to  increase,  becoalng 

larger  as  r becoaes  saaller.  As  r becomes  larger  than  r^  the  energy 

will  also  increase,  asyaptotically  approaching  the  sua  of  the  energies 
of  the  individual  atoas,  shown  as  the  zero  value  of  energy  in  the 
figure.  The  difference  between  this  zero  of  energy  and  the  energy  at 
the  average  intemuclear  separation,  r^,  is  the  dissociation  or  bonding 

energy,  E^.  An  aaount  of  energy  equivalent  to  E^  is  the  least  amount 

required  to  dissociate  the  molecule  into  its  constituent  atoas.  A 
similar  situation  occurs  when  a gas  approaches  a metal  surface;  the 
total  electronic  energy  will  vary,  the  occurrence  of  a minimum  isqplying 
that  the  gas  may  stick  to  the  metal  surface.  One  would  like  to  be 
able  to  calculate  whether  the  gas  does  stick  to  the  surface  and,  if  so, 
how  strongly  it  is  held  there.  The  characterization  of  these  problems 
requires  an  accurate  determination  of  equilibrium  separation  distances 
and  dissociation  qnergies.  The  importance  of  the  dipole  moment  is 
thought  to  be  in  its  relation  to  the  change  of  work  function  with 

crystal  synmetry  plane^. 

A method  which  we  have  explored  for  studying  these  various  aspects 

relating  to  chemisorption  is  the  X-o  method^,  a single-electron,  self- 
consistent-field  method.  The  tern  "single-electron"  means  that  the 
electrmis  are  each  considered  to  be  iMving  in  a field  resulting  from 
the  atomic  nuclei  and  an  average  charge  density  due  to  all  other 
electrons,  except  the  particular  one  in  question,  considered  together. 
This  is  in  distinction  to  a many-body  approach  where  the  interaction 
of  each  electron  with  each  other  electron  would  be  explicitly  taken 
into  account.  To  achieve  self-consistency,  one  starts  with  an  initial 
potential  for  the  aolecule,  V,  in  this  case  a siq;>erposition  of  atomic 


4.  L.B.  Shmidt,  in  Interactions  on  Metal  Surfaces.  Springer- 
Verlag,  Berlin,  lAV^,  p.  41. 

5.  J.C.  Slater,  The  Self-Consistent  Field  for  Molecules  and  Solids. 
Vol.  4 of  (hiantiai  ineory  of  Molecules  and  ^Hdfs,  Ncflraw-Hill', 
(Now  York),  1974. 
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potentials.  Pron  this  initial  potential  the  single-electron  eigenvalues, 
c^,  the  wave  function,  <|»,  and  the  charge  density,  p,  are  successively 

generated.  A new  potential  way  be  generated  from  the  charge  density. 

Hie  entire  process  is  iterated  until  the  new  potential  differs  insignifi- 
cantly from  the  ianediately  previous  potential.  When  this  occurs,  the 
calculation  is  said  to  be  self-consistent. 

For  a molecule  or  solid  the  Xa  theory  is  often  implonented  through 
the  multiple-scattering  formalism  which  adopts  a muffin-tin  (ffT) 

approximation^  to  the  potential.  This  muffin-tin  potential  consists 
of  a separate  spherically  symmetric  region  of  potential  centered  at 
each  atomic  nucleus,  regions  I and  II  in  Figure  2.  These  regions 
are  usually  chosen  to  be  touching.  The  entire  cluster  of  atomic 
Dotentials  is  surrounded  by  another  sphere,  outside  of  which  is  also  a 
region  of  spherically  averaged  potential,  region  IV.  This  sphere  is 
usually  chosen  just  large  enough  to  contain  the  atomic  regions.  In  the 
remaining  space,  which  is  neither  outside  the  large  outer  sphere  or 
inside  any  atomic  sphere,  region  III,  the  potential  is  taken  to  be  a 
constant.  Improved  agreement  with  experiment  has  sometimes  been  obtained 

7 

by  allowing  the  spheres  to  overlap  . 

This  type  of  potential  (MT),  while  quite  good  for  siany  purposes, 
is  too  severe  an  approximation  for  accurate  calculations  of  many 
properties.  To  remedy  this  situation  one  may  Incorporate  what  are 

g 

known  as  "non-muf fin-tin"  (hMT)  corrections.  It  would  be  costly  at 
present  to  calculate  the  "true"  Xn  potential  without  introducing 
restrictive  approximations.  However,  calculation  of  the  muffin-tin 
potential  and  the  first  few  terms  in  the  perturbation- like  expansion 
for  the  difference  between  the  "true"  Xa  potential  and  the  muffin-tin 
potential  is  not  exceedingly  costly.  Thus  one  is  able  to  calculate  a 
BK>re  accurate  approximation  to  the  Xa  potential.  These  NKT  corrections 
are  able  to  rewve  many  of  the  restrictions  of  the  isuffin-tin  potential 
while  not  increasing  computer-time  requirements  to  impractical  aisounts. 
The  NKT  corrections  consist  of  terms  «^ich  are  linear  (NMT-L)  in  Ap  and 
a term  which  is  quadratic  (NMT-Q)  in  Ap  i^ere  Ap  is  an  approxiution 
of  the  difference  between  the  x-o  and  x-a  muffin-tin  charge 
density.  The  energy  arising  from  these  NMT  terms,  < AE  > , is  given 


r. — ibfi,  pg.  loi. 

7.  K.H.  Johnson,  F.  Herman,  and  R.  Kjellander  in.  Electronic  Structure 
of  Polymers  and  Molecular  Crystals.  Plenum  (New  Vork),  1975. 

8.  J.B.  banese,  J.  biM.  t^kys.  61,  ^071  (1974). 
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«h«r«  V..  is  the  nucleer-electron  potential,  v is  the  electron-  electron 
no  c 

Coulo^  interaction  potential,  p is  the  charge  density,  and  a A preceding 

a term  indicates  a NKT  function. 

g 

In  previous  work  we  have  reported  the  results  of  calculations 
using  only  the  terns  of  A < > which  are  linear  in  Ap  in  a spin- 

independent  nodel.  The  present  calculations  extend  this  work  to  include 
the  tern  quadratic  in  Ap  utilizing  different  orbitals  for  different  spins 
(spin-polarized),  i.e.,  all  terns  in  (1)  were  calculated. 


III.  RESULTS 


A.  Excitation  and  Total  Energy  Calculations. 

Any  nethod  with  which  one  hopes  to  nodel  erosion  processes  in  gun 
tubes  nust  be  able  to  accurately  calculate  excitations  in  netals..gSince 
neay  radiationless  transitions  in  solids  have  quaslatonic  aspects*”,  we 
tested  the  ability  of  the  nethod  to  accurately  calculate  excitations  In 
netallic  atons  by  calculating  several  x-ray  satellite  lines  in  alunima. 
These  satellite  lines  occur  near  the  core  spectral  lines,  those  of  the 
inaemost  electrons,  as  the  outemost  electrons  undergo  traasitions. 

The  calculated  and  eq>erinental  results,  shown  in  Table  I,  indicate  that 
excitations  in  netallic  atons  nay  be  satisfactorily  studied  by  the  Xa 
nethod.  The  results  are  the  nore  striking  since  nultiply  ionised  atens 
are  being  described  by  a single  electron  theory.  The  value  of  a used  in 
each  case  is  that  which  vny  be  detemined  by  satisfying  the  virial 

theoran^^. 


Douglas  A.  Ringers,  "Electronic  Excitation  in  NO  and  8N  bjr  the 
Non-Nuffin-Tin-X-Alpha  Nethod",  ERL  Report  ItM,  US  An^  SalllatU 
Research  Laboratory,  Narch  1977.  (AD  fA037802) 

L.I.  Yin,  I.  Adler,  T.  Tsaag,  N.H.  Chen,  D.A.  Ringers  and  S. 

Crasenaan,  Phys.  Rev.  ^ 1070  (1974). 

Douglas  A.  Ringers,  "Xa  Calculation  of  Transition  Energies  ia 
Nultiply  Ionised  Atons",  IRL  NanorandiBi  Report  No.  2766,  June  1977. 

(AD  iA043091) 

j 1 - 

12  • -’t  ■ ’ 

V-”’  . - 

' ' -A 


Table  I.  Experiaentelly  Observed  X-Ray  Satellite  Lines  and  nieoretically 
Derived  Transition  Energies  in  Multiply  Ionized  Aluainua. 


n 

Experiment 

Theory 

Diagram 

Line 

Satellite 

Line 

Energy 

Difference 

(eV) 

Initial 
Vacancy 
Configumtion 
of  Satellite 

Pinal 
Vacancy 
Configuration 
of  Satellite 

Energy  Shift 
With  Respect 
To  Diagram 
Lines  (eV) 

•^2 

Ka' 

5.69 

l82p 

^P  2p,2p 

6.90 

•^2 

•^4 

11.85 

ls2p 

^D  2p,2p 

12.00 

19.71 

ls2p 

h 2p,2p 

19.63 

^2 

>^6 

23.56 

ls2p 

2p,3s 

23.15 

Aiamg  the  first  test  nolecules  studied  by  the  Xo  NMT  nethod  at  the 
Ballistic  Research  Laboratory  were  nitric  oxide  (NO) , carbon  nonoxide 
(CO),  and  sulfur  nitride  (SN).  The  open  shell  nolecule  nq 
has  been  the  subject  of  nuswrous  experinental  studies.  This  compound 
is  of  fundanental  inportance  in  ballistics  processes.  It  is  a product 
of  many  gas-phase  reactions  involving  propellants  and  explosives,  is 
corrosive  to  aetals,  and  is  a product  of  nany  cheailuainescent  reactions. 
Since  additionally  a great  aaount  of  experinental  infomation  is  avail- 
able on  NO  it  is  a desirable  molecule  upon  which  to  test  the  nethod. 

Much  of  the  sane  rationale  applies  to  carbon  nonoxide,  but  fron  a 
theoretical  or  nodelling  viei^int  it  differs  in  that  all  electronic 
occupied  orbitals  are  fully  ocoqpied  (closed  shell)  in  the  ground 
state.  Studies  were  also  perfomed  upon  the  nolecule  SN,  since  (SN) 

13  * 

is  an  explosive  , one-dinensional  conductor  and  is  sinilar  to  NO  in 

electronic  structure. 

i The  calculations  were  perfomed  in  a spin-polarized  node  (different 

I orbitals  for  different  spins)  for  the  MT  and  NNT-L  corrections.  Due 

to  their  snail  size  the  NMT-Q  corrections  were  calculated  in  a aonspin- 
polarized  nodel.  Partial  wave  expansions  up  to  i ■ 2 ("d  functions") 
were  used  on  all  centers.  The  a values  were  those  of  the  correspond- 
i ing  atom  (spin-polarized)  in  ths  case  of  open  shells)  in  the 

j atomic  regions  with  a suitable  average  elsewhere.  The  computational 

parameters  are  listed  in  Table  II. 


m Ford,  O.R.,  Phys.  Rev.  41,  577  (1952). 
13.  Shame,  J.,  private  conmunications. 
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Tabl*  II.  Calculatlonal  ParaMtws. 


Molacula 

HI 

“n 

“o 

“c 

“s 

“is 

NO 

2.17S8 

.74728 

.74118 

.74422 

CO 

2.17 

.74118 

.71881 

.74728 

SN 

2.8765 

.74952 

.72188 

.78274 

The  test  calculations  on  the  NO  nolecule  indicated  that  the  linear 
non>auff in-tin  terns,  E, , required  little  additional  computer  tine  beyond 
a nuffin-tin  calculation  and  amounted  to  1 percent  of  the  total  calcu- 
lated electronic  energy  of  the  molecules  (Table  III).  This  is  on  the 
order  of  the  bonding  energy  for  many  molecules  and  so  is  indeed  important. 
The  quadratic  term,  Eq,  took  considerably  longer  and  accounted  for  about 

0.1  to  0.2  percent  of  the  total  electronic  energy  (Table  III)  for  NO, 
an  amount  which  can  be  ignored  for  many,  but  perhaps  not  all,  excitation 
processes  as  calculated  by  total  energy  differences.  Nhen  studying 
excitation  processes,  one  is  interested  in  total  energy  differences  and 
the  difference  in  quadratic  texms  between  the  ground  state  (gs)  and 
any  excited  state,  i,  is,  for  NO, 

I Ep(gs)  - Ep(i)|  < .07  ryd  . 

If  we  consider  only  valence  states,  v, 

I E^Cgs)  - Eq(v)|  < .02  ryd  . 

Similar  results  were  observed  for  SN.  These  results  indicate  that  for 
most  excitations  quadratic  corrections  will  be  on  the  order  of  one 
eV  or  less. 

Since  the  quadratic  (NNT-Q)  term  was  the  most  computationally 
demanding,  an  effort  to  determine  if  there  existed  any  predictable 
dependence  upon  calculational  parameters  was  in  order.  Such  a depend- 
ency would  enable  the  NKT-Q  terms  to  be  estimated  for  a desired  set  of 
paraawters  from  only  a few  well-chosen  calculations. 

In  the  study  of  NO  the  outer  sphere  radius  was  increased  from  a 

value  corresponding  to  the  experimental  separation  of  2.1758  a.u.^^, 
where  the  outer  sphere  touches  the  atomic  spheres  and  is  just  large 

1^^  6.  Hertberg,  Molewlar  Spectra  and  Molecular  Structure,  I. 

Spectra  of  Diatomic  Molecules,  Van  Nostrand.  Princeton.  iSSO. 
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TABLE  in.  BINDING  ENERGY  CONTRIBUTIONS  IN  NO  (ryd.) 
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'•■S 


ii 


lOt  216.0862  3.0736  218.7'I62  39.7A  39.99 

Gno.  255.8908  2.9376  .3492  258.4842 


•nough  to  enclose  then  coapletely,  to  a value  .4  a.u.  larger.  The 
increase  was  in  incresients  of  .OS  a.u.  with  atoaic-sphere  radii  held 
fixed.  The  calculations  included  both  the  tenss  which  are  first  order 
in  dp  and  those  which  are  second-order.  The  dependency  of  the  quadratic 
tern  upon  outer  sphere  radius  is  approxinately  linear  as  shown  in 
Figure  3.  The  circles  represent  calculations  where  partial  waves  up 
to  t ■ 2 were  used  on  all  centers  and  the  squares  represent  calculations 
using  partial  waves  up  to  t ■ 3. 

For  a fixed  outer-sphere  radius,  the  ratio  of  atonic-sphere  radii 
was  allowed  to  vary  frosi  " .86  to  I^/Rq  " 1.17,  where  is  the 

radius  of  the  sphere  centered  on  the  nitrogen  nucleus  and  R^  is  the 

radius  of  the  sphere  centered  on  the  oxygen  nucleus.  The  quadratic 
correction  was  approxinately  linear  in  this  case  also  (Figure  4). 

The  total  energy,  however,  proved  to  have  a nonlinear  dependence 
iqpon  both  atonic-sphere  radius  (Figure  5)  and  outer- sphere  radius 
(Figure  6)  with  a definite  nlninua  in  both  cases.  The  ninimas  energy 
occurred  at  a value  of  Rji^/Rq  very  near  that  which  one  would  obtain 

using  the  respective  values  of  Slater  atonic  radii^^.  The  energy 

is  slight,  anounting  to  only  .15  rydberg  less  than  that  obtained 
for  an  outer-sphere  radius  corresponding  to  the  experinental  separation, 
and  occurring  about  .2  a.u.  larger  than  the  value  corresponding  to  the 
experinental  separation  distance.  The  values  of  atomic-  and  outer-sphere 
radii  corresponding  to  the  total  energy  isinina  were  used  in  ensuing 
calculations.  The  outer-sphere  radius  for  NO  of  2.3758,  iriiich  corresponded 
to  the  energy  nininun,  was  0.2  a.u.  larger  than  the  nlninun  radius 
required  to  enclose  the  atonic  region,  viz.,  2.1758  a.u.  This  larger 
value  is  narkedly  different  from  the  radius  usually  chosen,  i.e. , the 
nininun  radius  necessary  to  enclose  the  atonic  regions,  and  resulted 
in  i^>roved  results  for  dissociation  energy  and  electronic  excitation 
energies  in  NO. 

Using  these  selected  values  of  radii,  calculations  were  nade  of 
the  orbital  ionization  energies  in  NO.  The  ionization  energies  for 
each  occupied,  single-electron  orbital  were  calculated  as  total  energy 
differences.  In  Table  III  are  listed  the  nuffin-tin  energy,  linear 
NKT  correction,  quadratic  NNT  correction,  total  energy,  and  ionization 

potential.  For  comparison  experimentally  obtained  values^^  of  the 
ionization  potentials  are  also  listed. 


TJii  J.C.  Slater.  QuantSi  Theory  of  Molecules  and  Solids,  Vol.  2, 
McGraw-Hill,  New  York,  1965. 

16.  K.  Siegbahn,  et  al.,  ESCA- Applied  to  Free  Molecules.  North-Holland 
Amsterdam,  1971. 

* The  ordinate  is  taken  as  six-tenths  of  the  R^  to  R.  ration  for 
convenience. 
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Figure  3.  Quadratic  Energy  Correction  aa  a Function  of  Outer  ^ere 
Radius. 
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Ihese  results  are  plotted  in  Figure  7.  Note  that  the  1 ir 
level  Is  above  the  5 a level  of  corresponding  spin,  in  disagreenent 

with  reported  overlapping -sphere  calculations^^,  but  supported  by  other 

calculations^^.  In  the  nuffin-tin  calculation  without  MKT  corrections 
these  states  are  in  the  opposite  order.  For  this  simple  system 
omission  of  the  NMT  terms  yields  results  which  are  not  even  qualita- 
tively correct. 

Similar  calculations  of  the  ionization  energies  for  each  valence 
orbital  were  performed  on  SN  with  the  results  shown  in  Figure  8 and 
listed  in  Table  4.  In  SN  a situation  similar  to  that  in  NO  occurs: 
the  NMT  corrections  rearrange  the  ordering.  This  similarity  in  the 
ordering  of  the  ionization  levels  reflects  the  similar  ground  state 
electronic  structure  of  the  two  molecules.  The  ground  state  occupancy 
for  NO  is  a core  consisting  of  a closed  shell  (1  s shell)  plus  a 

valence  configuration  of  (3o)^  (4o)^  (So)^  (lir)^  (2h)^  while  that  of 
SN  consists  of  a somewhat  larger  core  of  closed  shells  plus  a valence 

configuration  of  (So)^  (6o)^  (7o)^  (2n)^  (Sir)^. 


Table 

IV. 

Binding  Energy  Contributions 

in  SN 

(ryd.) 

ORB. 

-^L 

-E.p 

Ip (calc. 

3w+ 

900.5439 

2.8089 

.4355 

902.9173 

.63 

2irt 

900.1338 

2.8355 

.4233 

902.5460 

1.00 

2w+ 

900.1980 

2.8409 

902.6156 

.93 

7of 

900.2938 

2.5936 

.4090 

902.4784 

1.07 

7er+ 

900.3719 

2.5385 

902.5014 

1.04 

6of 

899.8342 

2.7212 

.4115 

902.1439 

1.40 

6a4 

899.8998 

2.7189 

902.2072 

1.34 

Sof 

899.2885 

2.5285 

.3804 

901.4366 

2.11 

5o4 

899.3562 

2.5220 

901.4978 

2.05 

Gnd. 

901.2125 

2.7791 

.4461 

903.5455 

TTT- 

inr 

Salahub  and  k.lP 

Messmer,  J. 

Chen. 

Phys.  64,  2039, 

1976. 

18. 

P.B. 

Cade,  W.H.  Huo 

and  J.B.  Oreenshields , Atomic  Data  and  Miclear 

Data 

Tables,  1^,  1, 

1975. 
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Note  that  although  the  total  energy  in  the  ground  state  of  SN 
is  nearly  3.S  tiaes  that  of  NO,  differing  by  over  600  rydbergs,  the 
respective  linear  and  quadratic  non>euffin  tin  corrections  differ  by 
less  than  0.2  rydberg. 

A siailar  calculation  on  CO  produced  the  results  listed  in  Table 
V and  shown  in  Figure  9. 


Table  V.  NMT  Teres  and  Eigenvalues  for  CO. 


Energy 

in  ryd. 

State 

-ECMT) 

'q 

-Er 

-CjCexp.)^® 

or" 

lt3.S08 

2.^68 

.390 

185.^46 

39.69 

306 

C Is 

201.439 

2.541 

.267 

203.713 

21.66 

21.77 

3 0 

220.846 

2.152 

.249 

222.748 

2.62 

2.82 

4 0 

221.838 

2.383 

.301 

223,919 

1.45 

1.45 

1 t 

221.907 

2.511 

.274 

224.150 

1.22 

1.24 

S 0 

222.295 

2.415 

.339 

224.370 

1.00 

1.03 

G.S. 

223.217 

2.472 

.319 

225.371 

• 

’ 

The  two  excitations  listed  in  Table  6 show  that  although  the 
■uff in-tin  is  quite  good  for  certain  excitations,  it  cannot  be  relied 
upon  to  give  reasonable  values  for  all  excitations,  whereas  with  the 
NMT  corrections  all  calculated  excitations  were  quite  satisfactory. 


Table  VI.  Excitations  in  NO  (rydbergs). 


Excitation 

NT 

NKT 

Exp* 

(U)^ 

(3o)^ 

(2*)^  - (Iw)* 

(3o)^ 

(2»)^ 

.56 

.54 

.53 

(1*)^ 

(30)^ 

(2ir)'  - (lir)^ 

(3o)^ 

(2»)^ 

.34 

.65 

.63 

* Average  over  values  deduced  fro*  eapirical  curves  as  given  in 
reference  20. 


19.  Bagus,  P. , et  al..  Solid  State  Com.  20,  S,  1976. 

20.  LeFebure-Brion,  H.  and  C.M.  Noser,  J.~Qieai.  Phys.  M,  2951,  1966. 
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B.  Dipol*  Moamt  Calculations. 

Pravlous  calculations  of  tho  dipolo  Boaont  using  the  auffin-tin 
potential  without  the  NMT  corrections  have  been  soaewhat  disappoint- 

ing^^*^^.  These  calculations,  so  far  as  can  be  deterained,  were  all 
done  with  an  outer-sphere  radius  so  chosen  as  Just  to  encoapass  all 
atoaic  sjdieres.  Ihe  present  calculations  indicate  that  the  dipole 
aoaent  depends  strongly  upon  outer-  and  atoaic-sphere  radii,  as  one 
would  intuitively  expect.  The  dependence  of  the  dipole  aoaent  iqpon 
the  ratio  of  atoaic  radii  in  NO  is  shown  in  Figure  10.  When  the 
"energy  selected"  value  of  sphere  radii  is  used,  the  dipole  aoaent  for 
NO  was  calculated  to  be  .146  Debye  coapared  with  an  experiaental  value 

of  .157^^*^^  Debye.  This  value  coapares  favorably  with  the  result  of 
.147  Debye  obtained  in  very-large-basis-set  calculations  which  included 

25 

configuration  interaction  Siailar  procedures  for  SN  and  00,  Table 
VII,  also  were  in  agreeaent  with  experiaent. 


Table  VII. 

Dipole  Moaents  of  NO, 

SN,  and  CO. 

iteterlal 

Dipole  Moaent 

Slope  at  Equilibriua 

Calc. 

Exp. 

Calc.  Exp. 

rro 

.156 

m 1.17  (avg) 

SN 

1.834 

1.80 

CO 

.12 

.112 

The  dependence  of  the  dipole  aoaent  upon  intemuclear  separation 
was  calculated  for  several  distances  in  NO.  The  result  of  1.3  Debye/ 

26 

bohr.  Figure  11,  is  compared  with  a Hartree-Fock  calculation  (1.9 

27-31  ■ 

Debye/bohr)  as  well  as  with  an  experiaental  average  of  five 
values  (1.17  Debye/bohr  with  standard  deviation  of  3. OS  Debye/bohr). 

The  effect  of  the  dipole  aoaent  is  thought  to  be  iaqportant  in  eiqplaining 

32 

the  variation  in  work  function  with  crystal  sjnaMtry  plane 
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21.  Li,  £.Ii.,  Tnt.  JnlT  Quant.  Choa.  SIO,  193  (1976). 

22.  Woodruff,  S.B.,  and  Wolfsberg,  N.77.  Chea.  Phys.  S5,  36B7  (1976). 

23.  Burrus,  C.A.  and  Graybeal,  Phys.  Rev.  109,  1553  (l9?8). 

24.  Neuaann,  R.N. , Astrophys.  J.  161.  779  (i9'70) . 

25.  Billingsly  II,  F.P.,  J.  Chea.'TKys.  ^ 864  (1975). 

26.  ibid. 

27.  Schwin,  B.  and  Clough,  S.A. , J.  Chea.  PHy*.  38,  1855  (1963). 

28.  Jaaes,  T.C.,  J.  Chea.  Phys.  40.  762  (1904). 

29.  Schwin  B.  and  Ellis,  R.B.,  jTOiea.  ^ys.  45,  2528  (1966). 

30  Alaaichel,  C.,  J.  Phys.  (Paris)  ^ S4S  (1?%). 

31.  Chandraiah,  G.  and  Cho,  C.lf. , J.'Hblec.  Spectrosc.  134  (1973). 

32.  Schaidt,  S.D. , in  Interactions  on  Vatal  Surfaces,  Springer-Verlag, 
Heidelb^,  p.  84,  1976. 
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C.  IncTMsinf  tlw  Aeomcy  of  Total  Pnorgy  CalculatiOM. 

The  calculated  excitations  were  in  afreenont  with  expfrlMnt. 

These  excitations,  however,  are  calculated  as  differences  in  total 
energy  between  a ground  and  an  excited  state.  Although  this  energy 
difference  nay  be  sufficiently  accurate,  the  actual  energy  of  both 
states  My  be  in  error  by  a similar  saount.  For  excitations,  where 
there  would  be  an  error  cancellation,  this  would  be  of  little  cencefn. 

On  the  other  hand,  when  calculating  dissociation  energiea,  one 
needs  an  accurate  value  of  the  energy  of  the  ground  ftate  of  the 
nolecule.  Now  the  energy  difference  is  between  the  molecular  ground 
state  and  the  energy  of  the  individual  atoms  when  no  longer  found  in 
a nolecule,  i.e.,  at  infinite  separation.  Thus  there  is  no  longer  the 
possibility  of  an  error  cancellation  when  conparing  two  nolocular 
states  both  of  which  My  bo  in  error  by  a sinilar  anount.  The 
dissociation  energy  was  disappointing  even  with  the  calculated  energy- 

selected  radii  in  which  case  it  was  36%  of  the  experinental  value^’. 
(Nevertheless,  it  waf  a definite  inprovsMnt  over  the  value  of  only  12% 
of  experinent  which  was  obtained  with  an  outer-sphere  radius  corrofponding 
to  the  etpiilibriun  separation  distance.)  Increasing  by  one  the  nunher  of 
partial  waves  used  in  the  expansion  resulted  in  a different,  nore  negative 
total  electronic  oMrgy.  Varying  the  outer- sphere  radiue  to  doterniM  if 

the  nininun  in  energy  still  occurred  at  the  ssm  value  of  outer-sphere 
radius  showed  the  total  energy  decreasing  (Figure  12)  until  the  ^niUM 
outer-sphere  radius.  Just  large  enough  to  contain  all  atoeiic  spheres, 
was  used.  The  dissociation  energy  for  NO  was  now  72%  of  oxperMMnt, 
Inclusion  of  one  Mre  tern  in  the  partial  wave  eiq^sion  i^weved  this 
to  98%  of  the  ei^erinental  value  for  NO. 

The  saM  procedure  was  used  for  CO  with  the  results  as  shewn  in 
Table  VIII. 


Table  VIII.  Calculated  Bquilibriun  Separations  and  Dissociation 
Energies  for  NO  and  GQ. 


Conpound 

Dissociatioh  Eneriy 

tsv) 

BXP. 

NO 

2.175 

2.162 

6.S 

6.4 

GO 

2.142 

11.1 

11.5 

3Ti  G.  Hertberg,  op  ciV. 
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IV.  CONCLUSIONS 


Nuiy  calculated  valuaa  of  excitation  energies  and  ionisatloa 
potentials  utilising  the  nultiple-scattering  x-alpha  aethod  hetve  been 
in  good  agreenent  with  experinmts.  The  inclusion  of  the  MCT  correctiena 
which  have  the  ability  to  account*  to  a certain  extent*  for  aon-nuffla>tta 
portions  of  the  potential  hat  increased  the  usefulness  of  the  netbed. 

Ihe  non-nuffin>tin  CNKT)  corrections  have  predicted  binding  fbr  aolecules 
which  were  previously  unbound  in  calculations  using  the  nuffin-tin 
potential . 

The  nost  inporteat  aspect  of  the  present  work  it  that  a nethod  was 
developed  for  perfoming  accurate  total  energy  calculations  in  the 
x-alpha  nethod  without  the  usual  arbitrariness  in  gelection  of  pera> 
neters.  One  needs  tinply  to  use  relatively  high  t-values  (19  to  i*4  in 
NO,  for  exsiqple)*  choose  the  atonic  spheres  in  the  ratio  of  their  Slater 
atonic  radii*  and  use  an  outer  sphere  radius  equal  to  the  intemuclear 
separation. 

A systenatic  study  of  the  dependence  of  the  energy  \q>on  sphere 
radii  and  i-values  has  yielded  relationships  for  optinal  choice  of 
sphere  radii  should  the  use  of  large  1-values  be  undesirable.  This 
study  has  also  suggested  a nethed  for  detemining  dipela  noaents  of 
iaproved  accuracy. 

The  practical  i^lication  of  all  this  is  chat  Che  nethod  now 
appears  viable  for  ^e  study  of  chenisorption  processes  so  long  as 
one  includes  NNT  correction  terns.  A surface  could  be  nodelod  as  a 
cluster  of  netallic  atons,  as  schenatically  illustrated  in  Figure  IS* 
while  a gaseous  spocies  is  allowed  to  approach.  An  accurate  detomina- 
tion  of  the  potential  energy-curves  for  different  chenisorption  sites 
would  yield  the  forces  upon  the  gas.  For  different  gaseoua  species, 
atons  and  nolecules*  one  could  than  detemine  relative  nstal-gae 
adsorption  energies,  preferred  adsorption  sites*  and  in  the  case  of 
aolecules*  preferred  orientation.  An  accurate  calculation  of  dissociation 
energies  is  essential  to  the  study  of  chenisorption  and  the  chenical 
interactions  of  natal  surfaces  with  gaseous  species.  The  calculated 
values  of  dissociation  energy  for  NO  and  00  provide  an  inpetus  for 
applying  the  aethod  to  aetal-gas  systens  tdiere  the  ability  to  obtain 
accurate  values  of  these  properties  would  be  a significant  step  in  tile 
creation  of  a detailed  nicroscopic  theory  of  erosion. 
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